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Abstract 

A fast restoration of chiral symmetry in excited mesons is demonstrated. A minimal "realistic" 
chirally symmetric confining model is used, where the only interaction between quarks is the linear 
instantaneous Lorentz-vector confining potential. Chiral symmetry breaking is generated via the 
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, nonperturbative resummation of valence quarks self-energy loops and the meson bound states 

are obtained from the Bethe-Salpeter equation. The excited mesons fall into approximate chiral 
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multiplets and lie on the approximately linear radial and angular Regge trajectories, though a 
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. significant deviation from the linearity of the angular trajectory is observed. 
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There are certain phenomenological evidences that in highly excited hadrons, both in 
baryons QflQ and mesons |^,|5[ chiral and U(1)a symmetries are approximately restored, 
for a short overview see This "effective" restoration of chiral and U(1)a symmetries 
should not be confused with the chiral symmetry restoration at high temperatures and/or 
densities. What actually happens is that the excited hadrons gradually decouple from the 
quark condensates. Fundamentally it happens because in the high-lying hadrons the semi- 
classical regime is manifest and semiclassically quantum fluctuations of the quark fields are 
suppressed relative to the classical contributions which preserve both chiral and U(1)a sym- 
metries 0,0|. The microscopical reason is that in high-lying hadrons a typical momentum of 
valence quarks is large and hence they decouple from the quark condensate and consequently 
their Lorentz-scalar dynamical mass asymptotically vanishes . Restoration 

of chiral symmetry requires a decoupling of states from the Goldstone bosons 3, 0, 0, [l(3] 
which is indeed observed phenomenologically since the coupling constant for h* —> h + it 
decreases fast higher in the spectrum. 

At the moment there are two main paths to understand this phenomenon. In the first one 
one tries to connect the high-lying states to the short-rangepart of the two-point correlation 
function where the Operator Product Expansion is valid PjIuI- However, the OPE is an 
asymptotic expansion. Then, while the correct spectrum of QCD must be consistent with 
the OPE, there is an infinite amount of incorrect spectra that can also be in agreement 
with the OPE. Hence the results within the present approach crucially depend on additional 
assumptions 12. l^. 

In the second approach the authors try to understand this phenomenon within the micro- 
scopical models E.E, 10, wj]. There are also interesting attempts to formulate the problem 
on the lattice lfj, [l7|, though extraction of the high-lying states on the lattice is a task of 
future. 

In the absence of the controllable analytic solutions to QCD an insight into phenomenon 
can be achieved only through models. Clearly the model must be field-theoretical (in order 
to be able to exhibit the spontaneous breaking of chiral symmetry), chirally symmetric and 
contain confinement. In principle any possible gluonic interaction can contribute to chiral 
symmetry breaking and it is not known which specific interaction is the most important 
one in this respect. However, at the first stage it is reasonable to restrict oneselves to the 
simplest possible model that contains all three key elements. Such a model is known, it 
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FIG. 1: Graphical representation of the Schwinger-Dyson (A) and Bethe-Salpeter equations (B) 
in the ladder approximation. 

is a generalized Nambu and Jona-Lasinio model (GNJL) with the instantaneous Lorentz- 



vector confining kernel 
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This model is similar in spirit to the large N c 't Hooft 
model (QCD in 1+1 dimensions) 21] . In both models the only interaction between quarks 
is the instantaneous infinitely raising Lorentz-vector linear potential. Then chiral symmetry 
breaking is described by the standard summation of the valence quarks self-interaction loops 
in the rainbow approximation (the Schwinger-Dyson or gap equations), while mesons are 
obtained from the Bethe-Salpeter equation for the quark-antiquark bound states, see Fig. 
1. 

Conceptually the underlying physics is very clear in the 't Hooft model in the sense 
that once the proper gauge is chosen, the linear Lorentz-vector confining potential appears 
automatically as the Coulomb interaction in 1+1 dimensions. In 3+1 dimensions, once the 
Coulomb gauge is used for the gluonic field |22| , an almost linearly raising confinement 



potential has been obtained 



23, 2^ 



An obvious advantage of the GNJL model is that it can be applied in 3+1 dimensions 
to systems of arbitrary spin. In 1+1 dimensions there is no spin, the rotational motion of 
quarks is impossible, and the states are characterized by the only quantum number, which is 
the radial quantum number. Then it is known that the spectrum represents an alternating 
sequence of positive and negative parity states and chiral multiplets never emerge. This 
happens because in 1+1 dimensions the valence quarks can perform only an oscillatory 
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motion. In 3+1 dimension, on the contrary, the quarks can rotate and hence can always be 
ultrarelativistic and chiral multiplets should emerge naturally jsj. 

Restoration of chiral symmetry in excited heavy-light mesons has been studied with the 



quadratic confining potential 14j and it was also mentioned in a model with the instan- 
taneous potential of a more complicated form Q|. Here we report our results for excited 
light-light mesons with the linear potential. To our best knowledge this is the first explicit 
demonstration of the restoration of the chiral symmetry within the solvable "realistic" field- 
theoretical model. The results for the lowest mesons within the similar model have been 
previously reported in ref. j^J. 

The GNJL model is described by the Hamiltonian 19] 



H = j d 3 'x4)(x } t) (^—tf-x/ + mjilj(x,t) 

+ X -l d 3 xd 3 yJ^(x,t)K^(x-y)J b Mt), (1) 



with the quark current-current (J^(x,t) = ip(x, t)i^ip(x, t)) interaction parametrised by 
an instantaneous confining kernel K ab v [x — y) of a generic form. In this paper, we use the 
linear confining potential, 

Kf v {S-y) = g, g»o5 ab Vo(\x-y\), (2) 

and absorb the color Casimir factor into string tension, ^j-Vo(r) = ar. 

The Schwinger- Dyson equation for the self-energy operator E(p) = [A p — m] + ('jp)[B p —p\ 

is 

d A k - 1 



<E ® = I W¥ v{ v- kho s \k ,k) ^ 70 ' (3) 



where 



V(p) = - j d?xe^a\x\ = (4) 



so that the dressed Dirac operator becomes 

D(pa,p) = iS~ l {p ,p) = 70P0 - (lf)B p - A p , (5) 

where, due to the instantaneous nature of the interaction the time- component of the Dirac 
operator is not dressed. The Lorentz-scalar dynamical mass A p as well as the Lorentz-vector 
spatial part B p contain both the classical and quantum contributions, the latter coming from 
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loops 



B, 



A 



p ~ 




■V(p- k) siny2 fc , 
(pk)V(p- k) cos ip k , 



(6) 



(7) 



where tan ip p 




Solution of the Schwinger-Dyson equation (jSJ) with the linear potential is well-known, see 
e -g- jsiL an d the mass-gap equation has a nontrivial solution which breaks chiral symmetry, 
by generating a nontrivial dynamical mass function A p . This dynamical mass is a very fast 
decreasing function at larger momenta. Then the quark condensate is given as 



The homogeneous Bethe-Salpeter equation for the quark-antiquark bound state with 



where x(Pi M) is the mesonic Salpeter amplitude in the rest frame. Eq. (jOJ is written in the 
ladder approximation for the vertex which is consistent with the rainbow approximation for 
the quark mass operator and which is well justified in the large- Nc limit. 

The Salpeter amplitude can be decomposed into two components for mesons with J PC = 
(2n) h , (2n + 1) H , (2n + 1) ++ , (2n + 2) , or ++ and into four components for mesons 
with J PC = (2n + 1) or (2n + 2) ++ , respectively. Here J is the spin, P the parity and 
C the charge conjugation parity of the meson and n G No- In that way the Bethe-Salpeter 
equation becomes a system of coupled integral equations for the components which we solve 
by expanding them into a finite number of properly chosen basis functions. This leads to 
a matrix eigenvalue problem which can be solved by standard linear algebra methods. We 
vary the meson mass until one of the eigenstates is equal to one. 

In the gap as well as in the Bethe-Salpeter equations the infrared divergences are removed 
by introducing a finite "mass" into a confining potential, which is a standard trick. Then in 
the infrared limit ("mass" goes to 0) the quark propagator consists of a finite and diverged 
parts, while the mesons masses are finite. Recently it was demonstrated that also the masses 




(8) 



mass M in the rest frame, i.e. with the four momentum P M = (M, P = 0), is 




(9) 
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of quark-quark subsystems in the color-antitriplet state go to infinity in this limit and hence 

n 

are removed from the physical spectrum [28]. The results presented in the following were 
obtained by calculating the infrared limit numerically, i.e. the quoted meson masses were 
extrapolated to the infrared limit from a few points with a very small but finite mass of 
the infrared regulator. It turned out that in this region of the small mass of the infrared 
regulator the squares of the meson masses depend almost linearly on the mass of the infrared 
regulator making the extrapolation reliable. The presented results are accurate within the 
quoted digits at least for states with small J and for states with higher J but small n. At 
larger J for larger n numerical errors accumulate in the second digit after comma. 

By definition an effective chiral symmetry restoration means that (i) the states fall into 
approximate multiplets of SU (2) l x SU(2)r and the splittings within the multiplets ( AM = 
M_i_ — M-) vanish at n — > oo and/or J —>■ oo ; (ii) the splitting within the multiplet is much 

n n n 

smaller than between the two subsequent multiplets |4|, |5|, |6( . 

The condition (i) is very restrictive, because the structure of the chiral multiplets for the 

. nn 

J = and J > mesons is very different |4|, U- For the J > mesons chiral symmetry 
requires a doubling of states with some quantum numbers in contrast to the J = states. 
Given the complete set of standard quantum numbers /, J PC , the multiplets of SU(2) L x 
SU(2) R are 

J = 



(l/2,l/2) a : 1,(T + 
(1/2,1/2), : 1,0 ++ 

J = 2k, k=l,2,... 



-> 0,(T + 

-o,o- + , 



(10) 



< — ► o,r 

< — ► 0,J^ 



(0,0) : 0,J- 

(l/2,l/2) a : 1,J- + 

(1/2,1/2), : 1,J++ M 0,J 

(0,1)0(1,0) : l,J ++ < >1,J 



(11) 



J = 2k-l, k=l,2,. 



6 



(1/2, l/2) a 
(1/2,1/2), 
(0,1)0(1,0) 



(0,0) 



0, J+ 

1, J + 



+ 



< — ► 0, J 



< — ► 0, J~ 



1,J 



< — ► 0, J + ~ 

<— > 1, J ++ 



1,J 



(12) 



Note that within the present model the axial anomaly is absent. Even so there are no 
exact U(1)a multiplets, because this symmetry is broken not only by the anomaly, but 
also by the chiral condensate of the vacuum. Then the mechanism of the U(1)a symmetry 
breaking and restoration is exactly the same as of SU(2)l x SU(2)r. Hence the effective 
restoration of SU(2) L x SU(2) R would automatically imply restoration of U(1)a and of 
U(2) L x U(2) R and vice versa. An effective restoration of the U(1)a symmetry would mean 
an approximate degeneracy of the opposite spatial parity states with the same isospin from 
the distinct (1/2, l/2) a and (1/2, l/2) b multiplets of SU(2) L x SU(2) R . 

Note that within the present model there are no vacuum fermion loops. Then since 
the interaction between quarks is flavor-blind the states with the same J PC but different 
isospins from the distinct multiplets (1/2, l/2) and (1/2, l/2) 6 as well as the states with 
the same J PC but different isospins from (0, 0) and (0, 1) © (1, 0) representations are exactly 
degenerate. Hence it is enough to show a complete set of the isovector (or isoscalar) states. 

In Table 1 we present our results for the spectrum for the two-flavor {u and d) mesons 
in the chiral limit. Clearly the model should not be taken seriously for the low-lying states 
where other gluonic interactions as well as the 1/N C corrections should be important. The 
purpose of the study is, however, to demonstrate that a solvable field-theoretical model in 
3+1 dimensions does exhibit the effective restoration of the chiral symmetry at large radial 
excitations n and large spins. The excited mesons fall into approximate chiral and U(1)a 
multiplets and all conditions of the effective symmetry restorations are satisfied. We observe 
a very fast restoration of both SU(2)l x SU(2) r and U(1)a symmetries with increasing J 
and essentially more slow restoration with increasing of n. 

When the chiral symmetry breaking Lorentz-scalar dynamical mass of quarks is zero, then 
there are independent Bethe-Salpeter amplitudes just according to the chiral representations 
(10)-(12). A finite dynamical mass plays a role of the off-diagonal matrix element and mixes 
the otherwise independent chiral Bethe-Salpeter amplitudes for the states 1 , 2 ++ , 3 , .... 
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TABLE I: Masses of isovector mesons in units of yfa. 



chiral j-PC radial excitation n 

multiplet 1 2 3 4 5 6 

(1/2, l/2) _+ 0.00 2.93 4.35 5.49 6.46 7.31 8.09 

(1/2, l/2) 6 0++ 1.49 3.38 4.72 5.80 6.74 7.57 8.33 

(l/2,l/2) 1+- 2.68 4.03 5.15 6.14 7.01 7.80 8.53 

(1/2, l/2) b 1— 2.78 4.18 5.32 6.30 7.17 7.96 8.68 

(0,1) (1,0) 1— 1.55 3.28 4.56 5.64 6.57 7.40 8.16 

(0, 1) (1, 0) 1++ 2.20 3.73 4.95 5.98 6.88 7.69 8.43 

(1/2, l/2) 2~+ 3.89 4.98 5.94 6.80 7.59 8.31 8.99 

(1/2, l/2) b 2++ 3.91 5.02 6.00 6.88 7.67 8.41 9.09 

(0, 1) © (1, 0) 2++ 3.60 4.67 5.64 6.51 7.31 8.06 8.75 

(0, 1) © (1, 0) 2— 3.67 4.80 5.80 6.68 7.49 8.23 8.91 

(1/2, l/2) 3 + - 4.82 5.77 6.62 7.41 8.13 8.81 9.45 

(1/2, l/2) 6 3~ 4.82 5.78 6.65 7.44 8.17 8.86 9.50 

(0, 1) © (1, 0) 3~~ 4.68 5.63 6.48 7.26 7.99 8.67 9.30 

(0, 1) © (1, 0) 3++ 4.69 5.66 6.53 7.32 8.06 8.75 9.39 

(1/2, l/2) 4T+ 5.59 6.45 7.23 7.96 8.64 9.28 9.89 

(1/2, l/2) 6 4++ 5.59 6.45 7.24 7.97 8.66 9.30 9.92 

(0, 1) © (1, 0) 4++ 5.51 6.36 7.15 7.88 8.56 9.20 9.80 

(0, 1) © (1, 0) 4— 5.51 6.37 7.16 7.90 8.58 9.23 9.84 

(l/2,l/2) 5 +_ 6.27 7.05 7.78 8.47 9.11 9.72 10.3 

(1/2, l/2) b 5- 6.27 7.06 7.79 8.47 9.12 9.73 10.3 

(0, 1) © (1, 0) 5— 6.21 7.00 7.73 8.41 9.06 9.67 10.3 

(0, 1) © (1, 0) 5++ 6.21 7.00 7.73 8.42 9.07 9.68 10.3 

(1/2, l/2) 6 _+ 6.88 7.61 8.29 8.94 9.55 10.1 10.7 

(1/2, l/2) b 6++ 6.88 7.61 8.29 8.94 9.56 10.1 10.7 

(0, 1) © (1, 0) 6++ 6.83 7.57 8.25 8.90 9.51 10.1 10.7 

(0,1) ©(1,0) 6— 6.83 7.57 8.26 8.90 9.52 10.1 10.7 
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A key feature of this dynamical mass is that it is strongly momentum-dependent and vanishes 
very fast once the momentum is increased. When one increases excitation energy of a hadron, 
one also increases a typical momentum of valence quarks. Consequently, the chiral symmetry 
violating dynamical mass of quarks becomes small. Hence the mixing of the independent 
chiral Bethe-Salpeter amplitudes becomes small. A given state in the table is then assigned 
to the chiral representation according to the chiral Bethe-Salpeter amplitude that dominates 
in the given state. 

In Fig. 2 the rates of the symmetry restoration against the radial quantum number n 
and spin J are shown. It is seen that with the fixed J the splitting within the multiplets 
AM decreases asymptotically as 1/y/n, dictated by the asymptotic linearity of the radial 
Regge trajectories. This property is consistent with the dominance of the free quark loop 
logarithm at short distances. 

In Fig. 3 the angular and radial Regge trajectories are shown. Both kinds of trajectories 
exhibit deviations from the linear behavior. This fact is obviously related to the chiral 
symmetry breaking effects for lower mesons. Note, that the chiral symmetry requires a 
doubling of some of the radial and angular Regge trajectories for J = 1,2,... This is a 
highly nontrivial prediction of chiral symmetry. For example, some of the rho-mesons lie on 
the trajectory that is characterized by the chiral index (0,1) + (1,0), while the other fit the 
trajectory with the chiral index (1/2, 1/2);,. The intercepts of the asymptotic angular Regge 
trajectories for mesons in the given and different chiral representations coincide. Hence the 
asymptotic rate of the symmetry restoration with J is faster than 1 / \/~J . 

The numerical result for the quark condensate is (qq) = (— 0. 231 y/a) 3 , which agrees 
with the previous studies within the same model. If we fix the string tension from the 
phenomenological angular Regge trajectories, then yfo w 300 - 400 MeV and hence the 
quark condensate is between (—70 MeV) 3 and (—90 MeV) 3 which obviously underestimates 
the phenomenological value. Probably this indicates that other gluonic interactions could 
also contribute to chiral symmetry breaking. Notice that the string tension in Coulomb 
gauge can be larger than the asymptotic one. Lattice results suggest a value about twice 
the value obtained here j^, 0]. This would increase the value for the condensate but on 
the other hand lead to an unrealistically small pion decay constant 



31 



32]. 



In the limit n —>■ 00 and/or J —>■ 00 one observes a complete degeneracy of all multiplets, 
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FIG. 2: Mass splittings in units of \fo for isovector mesons of the chiral multiplets (1/2, l/2) a and 
(1/2, l/2)b (circles) and within the multiplet (0, 1) © (1, 0) (squares) against J for n = (top) and 



against n for J = and J = 1, respectively (bottom). The full line in the bottom plot is 0.7^/a/n. 



which means that the states fall into 



[(0,1/2)0(1/2,0)] x [(0,1/2)0(1/2, 



representation that combines all possible chiral representations for the systems of two mass- 
less quarks This means that in this limit the loop effects disappear completely and the 
system becomes classical jsl 7|. 
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FIG. 3: Angular (top) and radial (bottom) R^egge trajectories for isovector mesons with M 2 in 
units of a. Mesons of the chiral multiplet (1/2, l/2) are indicated by circles, of (1/2,1/2);, by 
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